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We study the low-temperature transport properties of clean one-dimensional spin-1/2 chains cou- 
pled to phonons. Due to the presence of approximate conservation laws, the heat current decays 
very slowly giving rise to an exponentially large heat conductivity, k ~ e T / ' T . As a result of an 
interplay of Umklapp scattering and spinon-phonon coupling, the characteristic energy scale T* 
turns out to be of order Qd/2, where Od is the Debye energy, rather than the magnetic exchange 
interaction J - in agreement with recent measurements in SrCuO compounds. A large magnetic 
field h strongly affects the heat transport by two distinct mechanisms. First, it induces a linear 
spinon-phonon coupling, which alters the nature of the T — » fixed point: the elementary excita- 
tions of the system are composite spinon-phonon objects. Second, the change of the magnetization 
and the corresponding change of the wave vector of the spinons strongly affects the way in which 
various Umklapp processes can relax the heat current, leading to a characteristic fractal-like spiky 
behavior of k(T, h) as a function of h. 



I. INTRODUCTION AND MAIN RESULTS 

One-dimensional (ID) magnetic systems exhibit a va- 
riety of interesting phenomena signifying their quantum 
many-body nature. They have been, therefore, the sub- 
ject of intense theoretical and experimental study. A 
great number of spin chain models has been proposed 
and investigated, with various interaction ranges, spin 
representations and anisotropies, as well as with cou- 
pling to other degrees of freedom. The nearest neigh- 
bor spin - 1/2 Heisenberg model, in particular, plays an 
important role being exactly integrableji allowing a de- 
tailed analysis of its thermodynamic properties. 2 These 
properties are found to be generic: they do not differ 
essentially from low energy properties of other (non inte- 
grable) short-range spin chains. On the other hand, the 
Heisenberg model does not provide a generic description 
of transport properties, fntegrability entails the existence 
of an infinite number of conservation laws, which in turn 
imply a dissipationless transport of the elementary ex- 
citations (spin-1/2 quantum solitons commonly named 
spinons) even at finite temperatures, T > . As a re- 
sult, measurable dc-transport coefficients such as electric 
and thermal conductivities are expected to be infinite. 
In realistic systems described by spin chains, e.g. ma- 
terials consisting of weakly coupled chains (or quasi-lD 
structures) of magnetically interacting ions, the conser- 
vation laws of the ideal model are partially violated, and 
the transport coefficients become finite. If the violation 
is soft one may expect an unusually high thermal con- 
ductivity in such systems attributed to spinon transport. 
Indeed, an experimental investigation of heat transport 
in materials of the (Sr,Ca,La)CuO series have found a 
considerable enhancement of the thermal conductivity 
in the direction parallel to the chains. 4 Other types of 
systems involving quasi-lD, antiferromagnetic spin-1/2 
chains such as organic compounds^ are in principle ac- 



cessible to such measurements, however at this point not 
much data is available^ Heat conductivity by magnetic 
excitations has also been observed, for example, in the 
spin-Pcicrls compound 7 CuGe03 or in ladder systems^ 

In a recent experimental study of SrCuO compounds, 9 
Sologubenko et at report and analyze specific heat and 
thermal conductivity data. The excess thermal conduc- 
tivity (k s ) along the chains direction (obtained after sub- 
tracting the phonon contribution) is identified as a con- 
tribution of the spinons, and its dependence on the tem- 
perature T is fitted to an empirical formula, designed 
to account for scattering of the spinons by localized de- 
fects as well as Umklapp processes. A phonon Umklapp 
mechanisms is suggested as an interpretation of the ex- 
ponential factor exp[T*/T] which describes k s (T) in the 
range of temperatures 50 K < T < 200 K, with a char- 
acteristic temperature scale T* ~ 200 K close to Orj/2, 
where ©d is the Debye temperature. However, this in- 
terpretation leaves open some questions: As the system 
is at zero magnetic field (corresponding, in fermionic lan- 
guage, to half filling) what inhibits the heat current relax- 
ation to proceed by low-energy spinon-spinon Umklapp 
scattering which is not exponentially suppressed? Why 
is not the relevant energy scale to the subleading relax- 
ation processes the exchange interaction J/ ks ~ 2500/4T 
rather than Orj/2 ~ 200 K? Another interesting observa- 
tion can be gleaned by conparing the plots of In k(T) vs 
1 /T of the pure phonon contribution and of the spinon 
contribution^ We note that the slope of the former is 
twice that of the latter. As pure phonon Umklapp pro- 
cesses are goverened by G, the fundamental reciprocal 
lattice momentum, we need to understand why it is G/2 
that dominates the spinon processes. 

We shall find that a rather subtle interplay of (approx- 
imate) conservation laws and quantum dynamics under- 
lies the experimentally observed heat conductivity, and a 
sophisticated hydrodynamic field theoretic approach is 



2 



necessary to fully account for it. The arguments are 
rather general^: When a system possesses some con- 
served quantities P these may "protect" the current J 
from degrading (this occurs when the cross-susceptibility 
Xjp 7^ 0) leading to a pure Drude peak and infinite 
d.c. conductivity. When the conservation of the pseudo- 
momenta P is softly violated they will, instead, lead to 
very long time tails in the decay of the current J, since 
states with a finite pseudo-momentum P typically carry 
also a finite current J (when xjp ^ 0). Hence the compo- 
nent of the current "parallel" to P, J^ P = (Xpj/Xpp)P 
with Xjp = X.j p pi will therefore decay exponentially 
slowly. The presence of such approximately conserved 
quantities leads then to a natural hydrodynamic descrip- 
tion of the system where a separation of fast and slowly 
decaying modes takes place and a consistent scheme of 
calculation of the slow mode conductivities can be carried 
out in terms of a matrix of decay rates of these modes. 

These features emerge when we study the transport 
properties of clean ID spin chains in the framework of a 
generic model which accounts for the coupling between 
spin excitations and lattice phonons. The model is ana- 
lyzed close to its fixed point Hamiltonian - the Luttinger 
liquid - and the thermal conductivity k(T) is evaluated at 
low but finite T. In this regime the transport coefficients 
are dominated by the slow relaxation of certain approx- 
imately conserved currents due to irrelevant corrections 
to the fixed point Hamiltonian. The most important cor- 
rection terms of this sort, which are relatively efficient in 
degrading the conserved currents of the integrable Lut- 
tinger model, are found to be associated with Umklapp 
scattering terms. These include pure spinon as well as 
spinon-phonon scattering processes. In the absence of a 
magnetic field, the latter class is shown to be dominant. 
Our detailed calculation of k(T) originating from Umk- 
lapp scattering of spinons by 3D phonons agrees with 
the experimental results of Ref. In particular, it ex- 
plains the origin of the exponential factor exp[T*/T] with 
7 • - (-);, 2. 

We then proceed to investigate the effects of a (large) 
external magnetic field h. The induced magnetization 
and the corresponding change of the wave vector of the 
spinons strongly modifies the way how various Umklapp 
processes relax the heat current. This leads to a fractal- 
like spiky behavior of k when plotted as a function of 
magnetization at fixed T, where the spikes occur at spe- 
cific commensurate values of the magnetization. Further- 
more, the magnetic field induces a linear spinon-phonon 
coupling tunable by h. The coupling alters the nature of 
the fixed point: the elementary excitations of the system 
are composite spinon-phonon objects. As a consequence 
of this mixing, the Umklapp processes are also modified 
and the relevant energy scale T* (again of the order of the 
minimum of Qu and J) depends smoothly on h. These 
effects are experimentally accessible in spin chains with 
relatively low magnetic exchange interaction J. 

The paper is organized as follows: in Sec. ^we derive 
the low energy model for the spinon system in the pres- 



ence of coupling to 3D phonons, and discuss the leading 
irrelevant corrections and their significance for transport. 
In Sec. IIIII we present the calculation of the conductiv- 
ity tensor by means of a memory matrix approach, and 
derive expressions for the thermal conductivity 
function of the temperature T. Sec. I1VI is devoted to 
the study of k(T) in a finite magnetic field h, where a 
linear coupling of spinons and ID phonons is accounted 
for within a Luttinger model. Our conclusions are sum- 
marized in Sec. [V] In Appendix ^ we emphasize that 
boundary conditions and finite size effects play an im- 
portant role in the presence of (approximate) conserva- 
tion laws and discuss what quantities are measured in a 
typical heat conduction experiment. Details of the cal- 
culation of the memory matrix elements are given in Ap- 
pendix [BJ For convenience, throughout the paper we 
adopt units where h = fiB — = L 



II. MODEL FOR THE WEAKLY COUPLED 
SPINON-PHONON SYSTEM 

We wish to compute the low-temperature thermal con- 
ductivity of a system consisting of a parallel array of long 
antiferromagnetic spin chains embedded in a 3D lattice, 
and interacting with the lattice phonons. Let us begin by 
describing the spin system. A typical spin chain model, 
with finite range interaction, is given by, 



H = 2 E ^ ( S ^ s i + S ^ S t) 

N N 

where S^, Sf are spin- 1/2 operators at lattice site i, 
and h is an external magnetic field applied along the z- 
direction. The coupling is antiferromagnetic: Jjj > 0, 
translational invariant: = J%—j and of finite range: 
Ji-j — for i — j > p. 

The low energy dynamics of this class of models is 
described by a Luttinger liquid Hamiltonian. A fully 
fledged derivation would proceed via repeated RG trans- 
formations and yield, in principle, the fixed point Hamil- 
tonian (the Luttinger liquid) as well as all irrelevant op- 
erators around it. Instead of following this route, we 
shall employ a short cut and proceed via the Jordan- 
Wigner transformation allowing a fcrmionization of the 
spin degrees of freedom and subsequent bosonization. To 
illustrate it we consider the XXZ model corresponding to 
the choice of range p = 1, but our conclusions are valid 
for a generic model. Begin by introducing a fermionic 
representation of the spin operators. In terms of spin- 
less fermionic operators tpi the spin operators can be ex- 
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pressed as 



11 



57 



ipi exp 



i-1 

1 

2 



Ipilpi 



The Hamiltonian H xxz is mapped onto a model of inter- 
acting spinlcss fermions in ID, 



nin i+ i + ft m (2) 



with J 2 determining the interaction strength and h play- 
ing the role of a chemical potential. The first term in 
H X xz corresponds to a kinetic energy term 



H k = -jy"^cos(fca)y4y>fc 



(where a is the lattice spacing) , whose low energy excita- 
tions (the spinons) are dominated by momenta k in the 
close vicinity of the two Fermi points ±fcp with 



k F = -(l + M) 



(3) 



where M = 2(S Z ) w h/(%J) is the magnetization of the 
spin chain (normalized to 1). 

In the low energy limit the field operator tpj is approx- 
imated as ip(x = ja) e fx iI>r(x) + e~ lkpX ipL(x) where 
the right and left moving fields ipn(x), ipL,{x) describe the 
low lying excitations, those near ±kp respectively. Using 
this expression in the Hamiltonian H xxz and keeping the 
leading terms only, one finds: 



YLy 



H 



LL 



where 



-i(Ja) I dx(ip R d x tpR - V4AV>i) 
+J Z J dx{p\ + p\+ W l ) , (4) 

Pr/l — ^r/l^R/l- ^he Luttinger Hamiltonian 
Hll thus obtained is conformally invariant. It is the 
fixed point of H xxz , but the coupling constant appearing 
in it are valid only to first order in J z / J and h/J. 

One may diagonalize -Hll by changing to bosonic vari- 
ables: the field 4>{x) and its conjugate H(x), satisfying 
[<f>(x),U(x')] — iS(x' — x). The fermion fields are then 
given bjii£ 



i 



R{L) 



Jl±1 



where 9 is defined by d x 8 = irll. The advantage of the 
bosonic representation is that it allows to diagonalize the 



interaction in H xxz by means of a Bogoliubov rotation 
leading to, 



in which (to leading order in \h\/J and \J Z \/J) 



IT 



a , K 



(5) 



(6) 



For arbitrary \J Z \ < J the Luttinger model Eq. (JSJ) 
still captures the low energy physics of the spin chain, 
however its derivation is more subtle. The exact Bethe 
Ansatz solution of the XXZ model yieldsii for h = (for 
h see Ref. 13]) 



K 



1 



2(1-^008-1 [i]) 



(7) 



In particular, in the physically interesting case of an 
isotropic Heisenberg antiferromagnet J = J z yielding 
K = 1/2. 

As noted above the transport properties of the XXZ 
model are not generic. The infinite number of conserved 
charges which assure its integrability also lead to a pure 
Drude peak and infinite d.c. -conductivity even at finite 
temperature. However the Luttinger liquid Eq. JSJ also 
describes (in the long wavelength limit) more complex 
spin chain structures as ladders, "zigzag" chains or in 
general chains with any finite range interaction (as long 
as no spin gap emerges) with the Jjj, Jfj dependence of 
the parameters v and K given by model-specific combi- 
nations of the coupling coefficients. The difference be- 
tween integrable and non-integrable spin chains in the 
low energy limit is captured by the structure of the irrel- 
evant operators around the Luttinger fixed point. In the 
latter case the irrelevant operators appear with generic 
coefficients. Consideration of the irrelevant operators 
is of crucial importance for transport properties. Non- 
integrable models are expected^ to have a finite heat 
conductivity at T > OMii^i^ To compute it, however, 
the fixed point Hamiltonian iJ^L is insufficient by itself: 
it is translationally invariant and integrable, and there- 
fore all the currents (e.g. spin current, heat current) de- 
scribed by it cannot degrade, leading to an infinite d.c. 
conductivity. One must add to the fixed point Hamil- 
tonian all irrelevant operators around it, and compute 
the conductivity from the resulting effective low-energy 
Hamiltonian. This implies, in passing, that the heat con- 
ductivity is a singular function of irrelevant perturba- 
tions, requiring us to recast perturbation theory in terms 
of a memory formalism (see below) . We also note that a 
number of recent studies 1819 found an infinite heat con- 
ductivity in generic non-integrable models. We believe 
that these claims are a consequence of either numerical 
problems^ or the neglect of certain classes of irrelevant 
perturbations^ 

We proceed to analyze the various irrelevant perturba- 
tions. In a generic model, all perturbations allowed by 



4 



symmetry are generated when high energy modes are in- 
tegrated out. The symmetries relevant for the following 
discussion are spin-rotation 1Z Z around the z-axis, dis- 
creet translations by a lattice spacing T a , inversion V and 
time-reversal T. In terms of the spinlcss fermions, 1Z Z 
guarantees charge conservation, T a leads to momentum 
conservation up to reciprocal lattice vectors G — 2n/a 
and the transformation rules under V and T are given 
by 

V : ip L -> ipR, 4>R -* ipL (8) 
T : if> L -> ip' R , ip R -> i -> -i. (9) 

The operators consistent with the symmetries above 
will be further divided into two classes, depending on 
their role in in transport phenomena. The first class, 
H u , consists of Umklapp operators H% m , describing pro- 
cesses where n spinons are moved from the right to the 
left Fermi point (and vice versa), possibly picking up 
m units of lattice momentum G. It is these processes 
that underlie the degrading of the currents leading to 
finite conductivities. The other class, H irr , contains low- 
energy processes where the number of spinons around 
each Fermi point remains conserved. This class includes 
corrections from band curvature, e.g. J ip R d 2 if)R, or from 



finite-range interactions. They do not affect the conduc- 
tivities directly. A formal way to distinguish between the 
two classes is as follows. Consider the two operators 

Jo = N R - N L , (10) 

where Nr and Nl are the total number of right and left 
moving spinons, respectively, and the spinon translation 
operator 

P Ts = j '(fe^-i^fi + ^WW' (11) 

These are among the infinite number of operators con- 
served by Hll, but play a special role in what follows. As 
we show in the next section, the conservation of certain 
linear combinations of Jo and Pt is minimally violated - 
in comparison with all other currents, the decay rates are 
exponentially small at low T. The class of operators Hi rr 
consists of terms in the Hamiltonian which conserve Nr, 
Nl and are invariant under (continuous) translations, 
hence commuting with both, [Hi rr , Jo] = [Hi rr , Pr s ] = 0. 
The first class, H u , includes Umklapp operators which 
do not commute at least with one of them. 



For even n the leading terms are of the form 

Jdx{e tAk - x j[^ R (x+3a)Mx+3a) + h.c.} = J dx(e iAk ™ V 2n <^> + h.c.) (12) 



ttU _ U 
nm "nm 



where 

Afe nm = n2kp — mG. (13) 

Note that Ak nm depends on the magnetization through kp and Eq. For odd n and vanishing magnetic field h, 
time reversal invariance does not allow terms of the form Ijl2(l , therefore leading perturbations are 

. 71+1 

HZn = 9nm / dx[e iAk - x (^ R (x)^ R (x) + ^ L {x)^ L {x)) n iP R (x + ja)iJj L (x + ] a) + h.c.] (14) 



i=i 



5 " m - / dx{e lAk ™ x e l2n ^ x) d x (f>(x) + h.c.) . 



(2tto) 



Thus far we discussed in detail the low energy descrip- 
tion of a single spin chain. We now turn to consider the 
complete system (as investigated experimentally e.g. by 
Sologubenko et al.)S consisting of an array of spin chains 
interacting with 3-dimcnsional acoustic phonons. 

The Hamiltonian of the system is 

H = H s + Rl D + H s , p , (15) 

where H s describes the spin array, H^ D is the phonon 
Hamiltonian, and H s p the interactions between spins and 



phonons. The spin array Hamiltonian is simply given by 
a sum of chains of the form we just discussed 

Hs = J2 H s ( 16 ) 

a 

where a labels the spin chains (parallel to the a;-axis) and 

s 12 LL ' irr ' Zjnm ^nra ' 

The Hamiltonian H^ D describes the system of three 
dimensional acoustic phonons to which the array of spin 
chains is coupled. In the following, we will consider 
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mainly phonons describing deformations of the lattice 
parallel to the chains which is chosen as the a;— direction. 
The dynamics of these deformations, to be denoted by q, 
is described by 



H 



:>,D 



d 3 x 
~27 



(17) 



where [i denote the x, y, z directions, P and q are (appro- 
priately normalized) canonical phonon momentum and 
coordinate operators and w M are the sound velocities 
with v x — v p and v y — v z = v± assuming a tetrago- 
nal symmetry. Acoustic phonons describing vibrations 
perpendicular to the chains are omitted here but 
can easily be included. As we will neglect in the fol- 
lowing the weak phonon induced interactions between 
different spin-chains, we will need only the propaga- 
tor of phonons along a single chain J d 2 k±G(ui, k) = 

Jd 2 k^_ l/{u 2 + „ ln[(£J 2 + V 2 k 2 )/Q 2 DA where 

0_d_l is the Debye frequency perpendicular to the chains. 
The space-time form of the resulting phonon propagator 
along a chain at T — is G p (t, x) ~ 1/ (x 2 + v 2 t 2 ). 

The phononic and spin degrees of freedom couple in 
a variety of ways which depend on the symmetries of 
the underlying lattice and on the strength of spin-orbit 
coupling. Assuming either inversion symmetry or weak 
spin-orbit coupling, the dominant coupling arises from 
the dependence of the exchange couplings on the dis- 



tance of atoms Ji 



J(Ri-Rj) m J + a(d x q)J' 



0[d 2 q, (d^q) 2 , (<9 M q^) 2 ]. In the presence of a magnetic 
field, other couplings can become important as will be 
discussed in section llVl In analogy to the classification of 
spinon-spinon interaction, we classify the various spinon- 
phonon interactions H S)P = Hf£ + £ nm H^~ p (again, 
generated by integrating out high energy modes) into 
Umklapp and non-Umklapp operators by their commu- 
tation relations with the operators Jo and the translation 
operator Pt of both spinons and phonons 

Jo = £ / dx [lp 1 Ra 1pRa ~ tpL a tpLa\ 

Pt = - f d 3 xPd x q + J2PTs (18) 

J a 

where P^ s is the translation operator of spinons [Eq. 
(|llfl ] on chain a. As we will show below, spinon-phonon 
couplings HYp" which commute with both Jo and Pr, for 
example J (d(f)) 2 d x q, will not be able to relax the heat 
current completely. More important are again Umklapp 
terms H„£~ p . For even n, leading contributions (in a 
given chain) are of the form 



H, 



U,s—p 



q U,P 

'hi m 



(27Td) 

while for odd n one obtains 



dx[e iAknmX e an *d a .q + h.c] (19) 



jjU,s-p 



q U,P 

(27ra) r 



dx[e 



e* n *(d x q){B x 4>)+h.a\ 



where (f> = (f> a , the phonon field q(R a ,x) is evaluated 
on the corresponding chain a, and Afc„ m is given by 
(113(1 . Processes involving multiple phonons are sublead- 
ing. We would like to emphasize again, that all the 
spinon-phonon couplings discussed above are irrelevant 
by power-counting. But we have to keep them if they are 
the dominating processes to relax the relevant approxi- 
mate conservation laws. Indeed, as will be shown in Sec. 
IIIII the spinon-spinon scattering terms Eq. (|12fl , (|14|) for 
which Afc nm is finite are exponentially suppressed with 
respect to the spinon-phonon terms in the physically rel- 
evant case, where the Debye temperature Qd is much 
smaller than J. 



III. THE THERMAL CONDUCTIVITY FOR 
ZERO MAGNETIC FIELD 



We turn now to the computation of the transport prop- 
erties of the spin chains, assuming the parameters of the 
system to be compatible with the SrCuO compounds 
studied in Ref. 9]. Since the exchange coupling J in 
these materials is extremely high (over 2000 K), magnetic 
field effects are negligible even at the strongest accessible 
fields (of order a few tens of Tesla). Hence, through- 
out this section we set h = 0. We note that in other 
compounds where J is much smaller (e.g., organic spin 
chains), interesting magnetic field effects should be ob- 
servable, as will be shown in Sec. IIVI We assume in 
addition the presence of inversion symmetry (or the ab- 
sence of spin-orbit coupling), so that the spinons couple 
to phonons only via H SyP = HfJ + J2 nm H nm~ p [ witn 
Hnm~ p given by Eqs. fllj, (gO)]' 

The transport properties of the spin chains at low 
temperature - like the charge transport in (quasi-) one- 
dimensional metalsi£*A£iSi - are governed by the approxi- 
mate conservation of certain quantities P ri m (to be called 
below "pseudo- momenta" ). The exponentially slow de- 
cay of a given P nm will lead to an exponentially large heat 
conductivity for low T . As already mentioned earlier, the 
reason is that states with a finite pseudo-momentum P nm 
typically carry also a finite heat current: the component 
of the heat current "parallel" to P nm will therefore decay 
exponentially slowly. 15 

We now proceed to identify the pseudo-momenta. 
Both Jo and Pt defined in l|18fl decay rather 
slowly as they commute with all non-Umklapp terms 
Hll, Hi rr , Hp D , Hl T p . More important are certain lin- 
ear combinations 



P, 



Akr, 



2n 



-Jo 



(21) 



(20) 



which we call "pseudo-momenta" (P n o is the usual mo- 
mentum operator). The pseudo- momentum P nm fur- 
ther commutes with all Umklapp terms with the quan- 
tum numbers n and m (and integer multiples kn, km), 



6 



namely: 



nrn 3 ±A t 
nm , 11 • 



3D 
P 
U 



H, 



H 



LL 
U.s- 



H£p] = 



J — u — [ r nr. 



T U,s-p] 



M -\- JT 

kn,km ' kn,kml 



For example, in the case of vanishing magnetic field con- 
sidered in this section, where kp = G/4 [see Eq. Q], 
P21 = Pt obviously commutes with all translationally 
invariant terms with Ak nm = and therefore with all 
possible low energy processes. 

The fixed point Hamiltonian Hll being conformally 
invariant possesses an infinite number of conservation 
laws. However, compared to our pseudo-momenta these 
modes decay much faster when generic perturbations are 
added since they do not commute with all the (low- 
energy) terms which we have collected in H irr and Ufl ■ 
They are therefore not important in the following discus- 
sion. The same argument applies to Hll + ^21 which is 
integrable or any other particular combination that hap- 
pens to possess conserved charges. 

Having established the presence of the slowly decaying 
modes we now turn to the question of how to calculate 
perturbatively transport coefficients in a situation dom- 
inated by a few of these modes. The method of choice is 
the memory matrix approachi^SiiSSiS^ which is based on 
the idea that while the conductivity is a highly singular 
function of the various perturbations this is not the case 
for the matrix of decay-rates of the slowest modes in the 
system. The memory matrix approach is formulated in 
a vector space of slowly decaying operators, spanned in 
our case by Pt, Jq and the heat current Jq as we want 
to calculate the heat conductivity. For convenience, we 
use instead of Pt, Jq and Jq the operators Jt,J s and 
Jq with: 



Jt = v 2 Pt and J s = vKJq, 
where J s is the spin current. In bosonized form 

Jt = —v 2 / d 3 x Pd x q + ^ / dxH a d x (j) a 
J s = vK^^ J dxTI a . 



(22) 

(23) 
(24) 



The heat current Jq = J d?xjQ X (along the chain di- 
rection x) is determined from the continuity equation 
d^jQfj, + dtH = where Ji is the energy density. For 
low temperatures it is sufficient to include only contribu- 
tions from the fixed point, J d 3 xTi ~ Hi ow with 



Hi ow — H, 



LL 



and one obtains 



Jq = - d 3 xv 2 p Pd x q-J2 



dxv 2 T\ a d x (f) a . (25) 



Adding further contributions e.g. from Umklapp terms 
or band curvature does not affect results to leading order. 
Note that the operators Jq and Jt are intimately related 
- in fact, they differ by the relative weight of the spinon 
and phonon degrees of freedom associated with the differ- 
ent velocities. However, there is a significant distinction 
between them: while Jt oc Pt remains conserved un- 
der all translational invariant corrections to Hi ow [even 
those which mix spinons and phonons, like f (d<fi) 2 d x q\, 
Jq does not remain so. 

We now set up the memory matrix formalism in the 
space spanned by the slow modes Jt, J s , Jq- To do so 
we follow Ref. (23 and introduce a scalar product (A\B) 
on the operators of the theory, 



{A{t)\B) = T f 1/T d\(A{t)^B{i\)) 
Jo 



(26) 



Then the dynamic correlation function of the operators 
A and B is 



Cab^) = / die™ (A(t)\B) 



A 



B 



iT 
to 



— / dte luJt ([A(t),B}) 



(A\B) 



(27) 
(28) 
(29) 



with the Liouville operator C defined by LA = [H, A] . 

In the space spanned by Jt, J s and Jq the matrix of 
conductivities is therefore given by 



TV 



(30) 



where V is the volume of the system and p, q are either 
of T, s and Q. The heat conductivity k is given by (c.f. 
Appendix 



(31) 



and a ss can be identified with the spin conductivity. The 
matrix of static susceptibilities can be written as 



Xpq 



1 

TV 



(Jp\Jq 



(32) 



As argued above, the matrix of conductivities a has no 
good perturbative expansion. We therefore express it in 
terms of a memory matrix M defined by 



a(tu,T) - x{T)(ti(w,T)-iux{T) 
and explicitly given as22, 



M pq {uj) 



1 

f 



d t Jp 



QCQ 



Q 



9 t Jq 



X{T). (33) 



(34) 



Note that in the literature^ the memory matrix is usu- 
ally defined as Mx _1 . The operator Q in Eq. (|34|l is the 
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projection operator on the space perpendicular to the 
slowly varying variables J p , 



(35) 



This separation between fast and slow modes underlies 
the perturbative expansion of M to which we now turn. 

The perturbative evaluation of M is greatly simpli- 
fied by the observation that since [Hi ow , Jk] = (for 
fc = s,T,Q), the operators dtJk are already linear in 
perturbations around the low energy Hamiltonian Hi ow . 
Hence, when these coupling terms are included to leading 
order in perturbation theory, one can set C = Ci ow with 
£>iow = [Hiow, ■] and Q = 1 in Eq. The expectation 

values in Eq. (|32() are also computed with respect to the 
low energy Hamiltonian Hi ow . Under these approxima- 
tions, the expression for the memory matrix Eq. (|34|) can 
be written as 



T 



nm 



M„ 



3-p) 



(36) 



where M nm and M nm . s -p are matrices in the space of 
the slow modes with matrix elements given by, 



M pq 



M pq 

nm,s—p 



(FP;F"Z-{FP;F"}°^ 



(37) 



(pP .pi \0 _ IpP .pi \0 



luJ 



Here F p = i[J p ,H u ], {F p ;F q )l is the retarded correla- 
tion function calculated with respect to Hi ow , and simi- 
larly for F p _ p = i[J p , H u ' s ~ p ] (the indices n, m have been 
omitted for brevity). Note that all the matrices M nm , 
M nrritS -p are symmetric. The static susceptibility ma- 
trix (for aT <§C Up) is given by 



(38) 










TTVT 2 


TTVT 2 


3 , 


3 2 


TTVT 


7TVT 


3 


3 



(where the matrix indices p, q take the values s, T, Q). 

We are mainly interested in the d.c. thermal conduc- 
tivity k{T) (Eq. Q2J at u = 0), which can be obtained 
from Eqs. (JSE) and (£23) through JSHJ in the limit oj -> 0. 
We find, 



k(T) 



TT 2 V 2 T 3 



(M-^tt + 2{M- 1 ) QT + {M- 1 )qq 



(39) 



with 



(M) Pq = Y,(.9n m i) 2 M p J(&k nm ,T) , (40) 



iitnl 



M pq (Ak nm ,T)^ limM" 



For conciseness we introduced the index I: I — denoting 
MJ^, and I = 1 denoting M p 'j ns _ p (and similarly for the 
coupling constants g ). Note that we have only retained 
contributions of Umklapp operators in M®Q . There are 
further contributions arising from H l s r p which turn out to 
be subleading for vanishing magnetic field and are there- 
fore omitted here. They are, however, important in the 
case of a finite magnetization, cf. section llVl 

The characteristic Luttinger liquid behavior of the 
spinon system is reflected by the functional dependence of 
M pq (Ak,T) on T and Afc. Approximate expressions for 
these functions can be obtained analytically in the high 
T limit, T > «|Afc|, or the low T limit, T < v p \Ak\, (see 
Appendix [B] for a detailed calculation) . We first note 
that the conservation law l|22[l implies a trivial relation- 
ship between F s and F T , and consequently for any T 
and Afc 



M*f(Afc,T) 
M„7(Afc,T) 
M« T (Afc,T) 



vAk 
~2nK 
v 2 (Ak) 2 

An 2 K 2 
vAk 



AC (Afc, T), 

M'f(Afc,T), 
M n t(Afc,T). 



(41) 



We therefore need to compute directly only three types 
of functions: M$, M® s and M% Q . In the high T limit 
we get 



AO Afc, T) 
M„t(Afc,T) 
M% Q (Ak,T) 



rp2(n 2 K+l)-3 

M% (Afc, T) Ak 
M s n f{Ak,T)T 2 . 



(42) 



More interesting is the low T limit T <C v p |Afc|, in which 
we find 



M-f(A*,r) 

M„t(Afc,T) 
M% Q (Ak,T) 



(2nKv) 2 M nl (Ak 1 T) 

AfcM^(Afc,T), 

(Afc) 2 M^(Afc,T) 



(43) 



where in the last line we have used the additional as- 
sumption (T/u|Afc|) <C (vp/v). The expressions for 
M„z(Afc,T) are the following: For even n 



M n0 (Afc,T)«e- 



and 



ir 2 T 2 {n 2 K/2)vT \ 2 



aAk 



2 K-2 

(44) 
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M nl (Ak,T) » A 



a 2 v 



(27ra) 2 "v 2 T(Afc) : 



2(n 2 _ff-l) 



~2v) 



2n 2 K 



a\Ak\ Y 



exp 



Up|AA;| 



2T 



Mm(Ak,T) » A 



a 2 T 



(27ra) 2n u2 w (Afc) 2 V t> / 



2v J 



exp 



Vp\Ak[ 
2T 



(45) 



(A is a numerical factor). For odd n, M n i are given by the above expressions multiplied by a factor ~ (Afc) 2 (for 
I = 0) or ~ (T/v) 2 (for I = 1), in particular 



Note that the above exponential factors are always dic- 
tated by the smallest of the velocities involved, and in 
our case v p = v m m = min{v,v p }. The physical origin of 
this behavior is that the minimal energy cost of a process 
involving a momentum transfer of Ak is associated with 
initial and final states of the elementary excitations with 
energy v m i n Ak/2 each. Since in the system of interest 
to us v p -C v, the exponential factor in l|44|) dramatically 
suppresses the pure spinon contribution to the sum in Eq. 
(I40|l (in particular, for |Afe| ~ l/a the exponent becomes 
~ — (J/T)). However, among the particular Umklapp 
scattering terms for which Ak = (and hence the 'high 
T" limit H42fl applies), the spinon-spinon process domi- 
nates as it contributes the leading power of T. 

We now focus our attention on the low T behavior of 
the thermal conductivity Eq. I|39|) . Eq. © implies that 
for h = 0, kp = Ti/2a = G/4. This is a particular, com- 
mensurate value of the filling 2kp/G, in which case the 
Umklapp term n — 2, m = 1 does not involve a momen- 
tum transfer, i.e. Afc 2 i = 0. Due to the exponential fac- 
tor in Eqs. (|44|l and (|45|l . the sum in Eq. i|40|) is strongly 
dominated by terms with a minimal Afc. In particular, 
the leading contribution to M ss and M®® is the single 
term n — 2, m = 1, I = corresponding to Afc 2 i = 0, 
where M%}(0,T) and Mg Q (0,T) are given by Eq. g2J. 
The other matrix elements vanish for Afc = 0. In fact, the 
vanishing of , M^i (for any p, q — s,T,Q) reflects 
the fact that Jt oc P 2 i commutes with all low energy 
terms in the Hamilitonian. Their leading contribution is 
therefore associated with the next smallest Afc, i.e. the 
term n = 1, m = and 1 = 1. These are given by the 
low T approximation (|46|l with n = 1, Afc = G/2. As a 
result 

C*" 1 ^ - M?T(G/2,T) (46) 

which is exponentially diverging. In contrast, (M^ 1 )qq 
and (M~ 1 )qt are inversely proportional to Mqq, associ- 
ated with the relatively fast short-time relaxation rate of 
the heat current Jq. As a consequence they depend al- 
gebraically on T and hence are exponentially suppressed 
compared to M^ T . Inserting into Eq. (|39|) this yields 



(47) 



with 



v p G 



(48) 



and kq depending on the parameters of the spinon- 
phonon system and the typical Umklapp scattering 
strength g 2 (« - <T 2 ). 

How does this compare to experiments? An exponen- 
tial behavior of the spin contribution to the heat con- 
ductivity has indeed been observed by Sologubenko et 
al? in SrCu2 and Sr 2 Cu03 for temperatures above 50 K 
(below which scattering from defects seems to become 
important). It was emphasized by the authors that T* 
is of the order 0z?/2, where Qd is the Debye temper- 
ature. A precise comparison to our result would re- 
quire a detailed knowledge of the phonon velocities in 
these systems. However, if we neglect for simplicity all 
anisotropies of the phonons, the Debye temperature is 
given by ®r> ~ v p (6n 2 /a 3 ) 1 / 3 w 0.6 v p G and therefore 
T* w 0.4 Q]j in very good agreement with the experi- 
mental observation. 



n(h = 0) re kq — 



/ rp x 


2(1-K) 






) exp 


T 



In we have defined k to be determined from 

an energy-current correlation function. However, in the 
presence of exact or approximate conservation laws and 
for finite systems it is far from obvious that this is the 
quantity measured in a typical heat transport experi- 
ment. For example, if spin is exactly conserved, bound- 
ary conditions will imply that no spin-current will flow 
through the surface of the sample and therefore the heat 
current in the experiment has to be calculated under the 
boundary condition of vanishing spin-current. As ex- 
plained in detail in Appendix lAl this implies that cfqq in 
Eq. (|3T1) has to be replaced by <jqq — 0q s /<7 ss . However, 
we consider in this paper a different limit, assuming that 
the sample is much longer than typical length scales on 
which e.g. the spin does decay. Under these assumptions, 
Eq. is indeed valid - see Appendix lAl for details. 

A note added on July 12, 2005: an error in the form 
of the phonon propagator lead to a wrong power-law pref- 
actor in Eq. (4-7) - the power 2(1 — K) should be replaced 
by —2K. See Appendix C for details. 
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IV. EFFECTS OF A FINITE MAGNETIC FIELD 

We now consider the effect of a finite magnetic field h 
on the heat conductivity of the spin chain. The field will 
have two main effects: the first effect will be to modify 
kp [see Eq. |J3J] and hence Ak nm leading to a fractal- 
like structure of the conductivity as a function of the 
magnetic field: as h is varied the system passes from 
incommensurable to commensurable values (for which 
Ak nm = for certain values of n, m) leading to a strong 
variation of the conductivity (see below). Clearly for this 
effect to be measurable, the spin-spin coupling J cannot 
be too large: to be observable with accessible fields, J 
needs to be of the order of a few tens degrees Kelvin. 
The second effect of a finite h is to induce linear phonon- 
spin coupling by the field, which alters the fixed point 
Hamilitonian of the system. Such a coupling is possible 
as the magnetic field breaks time reversal invariance T. 
(A similar coupling can arise as a consequence of spin- 
orbit interaction in crystals without inversion symmetry, 
even when h = 0). For finite h the linear coupling arises 
from terms of the form (d x q)SiSj ss (d x q)M5S z where 
M = 2(S Z ) is the magnetization and 6S Z = S z — (S z ). 

To analyze this case, we focus for simplicity on 
a strictly fD geometry considering only longitudinal 
phonons traveling along the chain direction. We note, 
however, that much of our forthcoming predictions are 
expected to be qualitatively applicable to spin chains em- 
bedded in higher dimensional systems as well. 

For ID phonons the free Hamiltonian Eq. I|17|l reduces 

to 



H„ 



dx 
2^ 



[(ttP) 2 + (d x q) 2 



(49) 



The normalization of P and q here is chosen differently 
than in Eq. I|17|l , so that their dimensions are the same as 
fl and 4>, respectively. It thus has the form of a Luttinger 
liquid with velocity v p and Luttinger parameter K = 1. 
To this we add a spinon-phonon coupling term of the 
form 



//, „ = -u I —n, oO,q . 



(50) 



At finite magnetic field, uq grows linearly with the mag- 
netization and as a consequence can be controlled. For 
h = and in the absence of inversion symmetry V similar 
terms which couple linearly to P rather than d x q arise 
from spin orbit coupling. They also give rise to a mixing 
of modes (the roles of P and d x q can be interchanged in 
the analysis below). 

A term of the form (|50|) leads to new eigenmodes of 
mixed spinon-phonon excitations. The Hamiltonian 



H = Hll + H„ + H s 



(51) 



is still scale invariant and is the fixed point of the coupled 
spinon-phonon system. We turn to diagonalize it. It is 



useful to define the free boson fields <j> = 4>/\K and 
II = yKU, in terms of which Eq. 151fl can be written as 



/,/ = / ^ v[{7T fi) 2 + (d x 4>) 2 ] 

Z7T 



(52) 



(u = uo\/K). We then diagonalize H* using the 
transformation, 25 



q 



C -Siv/Vp) 1 ' 2 
Siv/vp) 1 ' 2 C 



and similarly for the canonical momenta 
C 



p) {s(v p / v y/ 2 

where (assuming v > v p ) 



S{v p /vf/ 2 
C 



Hi 

n 2 



(53) 



(54) 



C 



s = 



1 
i 

71 



u 2 



1/2 



1/2 



(55) 



U 2 



[(v 2 - v 2 p f + u 2 vv p 



1/2 



The transformation l|53l) . 154fl is symplectic in order to 
preserves the canonical commutators, [(f> u (x),H u i (x')] = 
i8 vv i8{x — x'). The resulting Hamiltonian H* takes the 
diagonal form 



h* = [yJ2 v » U^u^ + ^id^) 



where 

1 



— = — r 2 

K 2 ■>.. l v 



(56) 



U 2 



viKi 



v. 



2v p 

v 2 K 2 = v p . (57) 



The strongest mixing (for a given u) occurs when v ~ 
Vp, which is realized in spin chain systems provided the 
exchange interaction J is not too large, and comparable 
to Qd- We then define v = (v + v p )/2, Sv = v — v p , where 
Sv <C v. Assuming in addition u <C v, we obtain 



Ki 
where 



v[l + A(h)}, v 2 w v[l - A(h)} , (58) 



1 



2v 



A(h) 



2v 



[4(Sv) 



,211/2 



4v 



(59) 



which depends on the magnetic field h via the coupling 
u. 
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We proceed to evaluate the heat transport coefficient 
following the memory matrix method described in sec- 
tion lllll The linear coupling modifies the low-energy heat 
current Jq , which is now defined with respect to the ele- 
mentary degrees of freedom of H* rather than Hll + H p . 
In terms of the two eigenmodes of H*, Jq is given by 



Jq = - / dx vlH v d x <fr u 



(60) 



In analogy with Eq. (|23|) . Jt (now defined as Jt — v 2 Pt) 
is given by 



Jt = —v 2 / dx ^2 Tt u d x (f>v ■ 



(61) 



i/=l,2 



The derivations of the memory matrix and the static sus- 
ceptibility proceed in the same way as in Sec. IIIII In 
particular, the memory matrix M is generally given by 
Eq. igSJ and by gUl in the oj — > limit. The static 
susceptibility matrix \ is given by an expression nearly 
identical to H38[). except that in xtt, Xtq and xqq the 
velocity v is replaced by 2v (to leading order in A(h)). 
The most prominent modification compared to the h = 
case is indicated in the spinon-spinon contribution to M 
(i.e., the I = terms): as a result of the composite na- 
ture of the elementary excitations, it is no longer a pure 
spin contribution, and is dominated at low T by the same 
exponential factor as the I = 1 terms. Since the latter 
involve a higher power of T, M£ 9 (Afc„ m , T) = M^ ac- 
tually dominates for any Ak nm . We therefore neglect the 
/ = 1 terms, and Eq. (|40|) becomes 



M(T) = J2(9n m ) 2 M n (Ak nm ,T) 



(62) 



Similarly to Eq. (|4*T]l . we find 



M£ (Ak,T) 



v 2 Ak 
2nKv 



M-(Afc,T) 



M?/{Ak,T) = ^^M%(Ak,T), (63) 
M Q , T (Ak,T) = -t^ M % s (Ak,T). 

nl \ 'J 2nRv nl V . J 

The calculation of the functions MP q (Ak, T) (for p, q de- 
noting either of s, Q) is essentially the same as in the 
h = case (see Appendix [5J). In the high T limit 
(T > ui|Afc|) we get 

M^(Ak, T) ~ T ^ n +M-3 
M® s (Ak, T) ~ M s n s {Ak,T)Ak (64) 
M% Q {Ak,T) ~ M° s (Ak,T)T 2 , 



where 



n^XiC 2 , l3 n = n 2 KK 2 S 2 (v/v p ) (65) 



and C, S 1 are defined in Eq. Q55J1. The low T limit cor- 
responding to T < w 2 |Afc| and T/{v x \Ak\) < 5 < 1 
yields 

M s n s {Ak,T) = {2nKv) 2 M n (Ak,T) , 
M® 8 (Ak,T) ~ AkM n (Ak,T) , (66) 
M% Q (Ak,T) ~ (Afc) 2 M„(Afc,T) 



with 



M n (Ak,T) 



Aa 2 



{2na) 2n vT{ Ak) 2 [A(h)] a " (2/3 n ) 2 ^ 



aT\ 2an (a\Ak\\ 2pn 



^ 2 |Afc| 
2T 



(67) 



Here A is a numerical factor, A(/i) is defined in Eq. H59|) 
and we have used the fact that «2 = min{ui,W2}. Note 
that since time-reversal symmetry is broken at finite h, 
Eq. i|12|) is the leading Umklapp term for arbitrary n, 
and hence Eq. (|67|l hold for both odd and even n. 

Inserting Eqs. i|63[l through (|67|l into l|62|) we get 
the dominant contributions to most elements of the ma- 
trix M(T). An exception is Mqq, which includes addi- 
tional corrections neglected in the above approximations: 
these are associated with irrelevant perturbations such as 
J (d(f>) 2 d x q, which do not commute with the heat current 
Jq . These lead to contributions to Mqq which are power 
law in T . In the case of a finite magnetization, when all 
the Umklapp terms are either suppressed exponentially 
[Eq. (|67() ] or by large powers of T (|42() . these corrections 
cannot be neglected. Hence, at low T one always gets 
Mqq 3> Mtt,Mtq and k is dominated by M^ T . Us- 
ing that Mqq is exponentially larger than the smallest 
eigenvalue of the matrix M, we obtain with exponential 
precision M^t w M SS /(M 8S M TT - M 2 T ) and therefore 
for the thermal conductivity k 



, x 4ir 2 v 2 T 3 . 1( , 
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V 



(68) 



with 



V ~ \ X X! M nm M n , m .{nAk n , m , - n'Ak nm f (69) 



nm n m 



and M nrra as an abbreviation for (gK m ) 2 M n (Ak nm ,T). 

How does the conductivity /c depend on T and /i? 
As discussed above, a magnetic field h leads to a linear 
spinon-phonon coupling and therefore tunes the param- 
eters A(h) and v-i in Eq. (|67|l . However, for large fields, 
of the order /i ~ J ', another effect is even more impor- 
tant: the finite magnetization M induces a shift of the 
Fermi momentum fc_F according to Eq. (PJ). The filling 



2k F /G 



_ l+M 



is set to an arbitrary, generally irrational 



value, and can be tuned continuously by varying h. Upon 
tuning kp, the characteristic momentum transfer Ak nm 
(|13|l associated with an Umklapp process H^ m is mod- 
ified accordingly. As Umklapp processes at low T are 
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suppressed exponentially by e " 2Afe ™ 



,/(2T) 



a change in 



1/51/4 1/3 2/5 1/2 3/5 2/3 3/4 



Ak nm modifies exponentially the contribution of H% m to 
the various relaxation rates. 

We first analyze n for low T and an almost commen- 
surate magnetization with M « 2^ — 1 or k F w ^f 1 ^ 

o n r 2 n 

where mo and no are (small) integer numbers. In this sit- 
uation, H^ oma is the strongest Umklapp term as, accord- 
ing to l|13[l . Afc„ omo ss 0. However, due to the approx- 
imately conservation of the pseudo-momentum P nomo , 
this Umklapp cannot relax the heat current alone, i.e. 
T> in Eq. I|t)8|) vanishes if only no and mo are included in 
the sum. Hence k is determined by the second strongest 
Umklapp H%i with the smallest possible (but finite) 

™0 m 

momentum transfer, Ak n > m > = ±G/no (the correspond- 
ing values of n' ,m' depend strongly on no, mo). Wc 
therefore obtain for the heat conductivity close to com- 
mensurability, 



1 o exp 



v 2 G 



2n T 



(70) 



where a n is defined in Eq. I|t)5|) . The expression for n com 
is valid as long as kp is sufficiently close to the com- 
mensurate value so that Ak„ orno <C T/vi (in which case 
the "high T" approximation holds for M™(Ak nomo ,T)). 
Note that the conductivity n is largest close to a com- 
mensurate magnetization M where no is small in ap- 
parent contradiction to the expectation that Umklapp 
is most efficient for commensurate fillings. The reason 
is simple^: while the strongest Umklapp is enhanced for 
commensurate M, the second-strongest which determines 
the size of k is suppressed. 

How large is the conductivity k for a typical incom- 
mensurate magnetization or for temperatures where the 
asymptotic behavior l|70|) is not yet reached? Umklapp 
processes H^ m with small n are suppressed at low T be- 
cause Ak nm is large while contributions with large n are 
suppressed by algebraic prefactors with large exponents 
a n and f3 n [see Eqs. (|64H67fl ]. We therefore estimate k(T) 
using Eq. I|t)8|) as follows: first, the summation over to is 
performed for a given n, noting that the dominant term 
(which maximizes M nm ) is to = mo, where mo/n is the 
closest rational approximation of IkpjG. The second- 
strongest Umklapp (for this n) is therefore characterized 
by a momentum transfer Ak ss G/N where N = na 
(with a of order unity) . We then evaluate the remaining 
sum over n in a saddle point approximation. This yields 



Kt yP (T) ~ exp 



2/3' 



(71) 



where f * = C[ln(v 2 G/T)] 1 / 2 v 2 G and C is a constant of 
order unity. 

In Figs. and [21 we show schematically the dependence 
of \d[k] on the magnetization. For a precise quantitative 
prediction of k(M, T) the RG flow of all Umklapp terms 
has to be calculated. While this is in principle possible 
e.g. for a highly anisotropic spin chain (small J z ) this is 
beyond the scope of this paper. To obtain the schematic 
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FIG. 1: Schematic plot of the dependence of the logarithm of 
the heat conductivity k(M, T) on magnetization M for various 
temperatures, k is largest for almost commensurate magneti- 
zations and low temperatures. Parameters used for this plot: 
1, G = 2tt, v = 1, v p = 0.8, K = 0.7, u = 0.1M. 



9nm 



picture of Figs. Q] and [3 we have set all g^ m to unity and 
used the asymptotic expressions H64H67H for the memory 
matrix in (|68|) . As expected, the M dependence of k(M) 
becomes more and more "spiky" towards lower temper- 
atures with maxima close to commensurate magnetiza- 
tions, e.g. for M = 0, |, ^, |, ... Unfortunately, the mul- 
tiple peaks in k(M) occur at extremely high values of n 
where in most experimental systems the heat transport 
will be dominated by impurities or sample boundaries. 
Therefore this effect, while being an amusing theoretical 
prediction, is difficult to be observed experimentally. At 
high temperatures, however, a pronounced minimum in 
k should be experimentally observable (see Fig. |5J in a 
regime where inelastic scattering still dominates trans- 
port. The precise position of this minimum depends on 
the temperature and on details of the system under con- 
sideration. We also would like to point out that a pic- 
ture similar to Figs. ^ and |2 would emerge even in the 
absence of spinon-phonon coupling, only the relevant en- 
ergy scales would be different and set by J rather than 
the minimum of Qd and J. 

Note that the position of the maxima in Fig. ^ seems 
to be shifted away from the commensurate positions. 
Consider for example the situation when M is close to 
§, M = \ + SM. As |Afc„ m | = G|n±±^ - m\ ac- 
cording to Eqs. © and l|13|) . the dominating scattering 
process is H^ 2 with Ak% 2 — ^-SM being very small. 
Therefore the heat transport is dominated by the relax- 
ation of P32 by the second strongest Umklapp H^i with 
|Ajfe u |=G(|-|5M) andln«;« ^(|-±<5M)for small 
SM and low T. The position of the local maximum to 
the left of M = 1/3 is determined by the competition 
with H% x with |Afei| = G{\ + SM). While the expo- 
nential factors in Eqn. I|67|l favor for SM < as 
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10 k ... i l u I ^-L^ I i_U I u I I . ^ 

0.2 0.4 0.6 0.8 

M 

FIG. 2: Schematic plot of the magnetization dependence of 
the heat conductivity k for various temperatures. Tempera- 
tures are higher than shown in Fig. otherwise parameters 
are identical. 

|Afe 2 i| - |Afcn| = ^j-SM < 0, the higher order process 
H21 is suppressed by algebraic prefactors and will there- 
fore only prevail for sufficient large 8M < or sufficient 
low T. 

In comparison to Ref. it is interesting to note that 
while the electrical (or spin-) conductivity in Id systems 
is suppressed for exactly commensurate fillings, this is 
not the case for k. The reason is that the overlap of the 
relevant pseudo-momentum and the heat current is finite 
at commensurate fillings, while the corresponding overlap 
of P-nm and spin or charge currents vanishes with expo- 
nential precision^ - therefore also the Wiedemann- Franz 
law will be violated exponentially for commensurate fill- 
ings. 



V. CONCLUSIONS 

In this paper we have studied the thermal conduc- 
tivity k of clean spin-chains coupled to phonons. The 
heat transport in the absence of defects is strongly in- 
fluenced by approximately conserved pseudo-momenta. 
Due to their presence, low energy processes cannot relax 
the heat current and therefore K is exponentially large at 
low T. The exponent is determined by the slowest mode 
in the system, i.e. in most materials by the phonon ve- 
locity. In semi-quantitative agreement with experiments 9 
we find k ~ exp[6_D/(2cT)] where Op is the Debye tem- 
perature and c = ^pg- a number close to 1 depending 
on the anisotropies of the phonon system. It is instruc- 
tive to compare this to the pure phonon Umklapp which 
leads to a relaxation of the phonon contribution to the 
heat current characterized by twice as large an exponent, 
k ~ exp[— v p G/(2T)]. The reason is that in the mixed 



spinon-phonon process considered by us the relevant mo- 
mentum transfer is G — 2kp — G/2 rather than G. In- 
deed, by comparing the plots of In k(T) vs. 1/T of the 
phonon and the spinon contribution in Sr^CuOj— one 
can identify slopes differing by this factor of 2. 

In the presence of a magnetic field, spinon and phonon 
modes start to mix. Furthermore, the Umklapp processes 
depend exponentially on the magnetization. This leads 
to a fractal-like spiky dependence of k on magnetization. 
Surprisingly, k is largest for commensurate magnetiza- 
tions. This is again a consequence of the approximately 
conserved pseudo-momenta: it is not the strongest, but 
the second strongest Umklapp process which determines 
the thermal transport. 
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APPENDIX A: CONDUCTIVITY, 
APPROXIMATE CONSERVATION LAWS AND 
BOUNDARY CONDITIONS 

In this appendix we discuss on general grounds, how 
boundary conditions and approximate conservation laws 
influence transport measurements. Let us assume that 
in an experiment the transport of a conserved charge qi 
(e.g. the energy density) is measured using a 4-point 
probe. We want do investigate the influence of a sec- 
ond approximately conserved charge q 2 (e.g. the spin- 
density) on the experiment. The arguments given below 
can easily be generalized to include further exact or ap- 
proximate conserved charges. The relevant continuity 
equations read 

d t qi+d x]1 = (Al) 
d t q2 + d x j 2 = -qz/T2, (A2) 

where r 2 describes phenomenologically the (slow) relax- 
ation of q 2 (assuming that q 2 = in equilibrium by def- 
inition). To set up a hydrodynamic description of the 
measurement, we assume that currents are driven by the 
force Fi (e.g. d x T/T) and by gradients of q 2 

ji = <r 11 F 1 - D 12 d x q 2 (A3) 
32 = 021F1 - D 22 d x q 2 , (A4) 

where the conductivity <j\ 2 describes e.g. the spin-analog 
of thermopower. In steady state, dtqi = 0, one obtains 
easily from Eqs. HA1HA4|) 

h = 3\°2\I<J\\ - D 22 d x q 2 (A5) 
D 22 d 2 x q 2 = q 2 /T 2 (A6) 
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with ji(x) = const, and D22 — -D22 — -Di2°"2i/< 7 n- These 
equations have to be solved with the appropriate bound- 
ary conditions. For our example the experimentally rel- 
evant boundary conditions are j2(±L/2) = (no spin- 
current flowing out of the sample) where L is the size 
of the sample. From l|A6p one therefore obtains for 
-L/2 <x< L/2 



<?2 = q° I exp 



x + L/2 
h 



exp 



x-L/2 
h 



(A7) 



where q® is determined from the boundary conditions and 
1-2 = V D22T is the (spin-) diffusion length. 

Obviously, one has to distinguish two different situa- 
tions. If L ^> I2 both q2 and d x q2 vanishes exponentially 
in the sample. Therefore according to (|A3|> the (heat) 
transport of qx inside the sample is determined by axx- 



jx = cfuFi for L > l 2 . 



(A8) 



This is the situation considered in this paper (formally, 
within our Hamiltonian l 2 — 00 as S z is conserved, but 
in a real material S z will decay e.g. due to spin-orbit 
scattering from impurities). Note that Eq. (|A4|) implies 
that the (spin-) current j'2 inside the sample is finite for 
o~i2 7^ 0, i.e. for h 7^ 0. This current, however, decays 
close to the sample boundaries in a width of order I2. 

The situation is different in the second case when L <C 
I2 or if q2 is exactly conserved. Then j% vanishes inside 
the sample due to the boundary conditions and plugging 
[H into (IA3I) one obtains 



n 



-D12C21 



D 



22 



Fx forL<Z 2 . 



(A9) 



order in (v p /v) < 1, we find 



™ M _ -Wnml 2nKv 

^nmA 1 ) - 



(2na) n 



(Bl) 



F nml (t) 



where for abbreviation we have omitted the indices n, m 
from Ak nm and introduced the definition 



{2ira) n 

x I dx{e l ^ kx e l2n ^ x ^d x (t>bi{x) - h.c.) , 



bi(x) = (d x q) 



(B2) 



We then find that the retarded correlation function 
( F nmV F nml)l at frequency w is given by2i 



( F nmV Fnml)u — (^3) 

/oc />oo 
dx \ dt e l M- A M lm{C pq (x,t)} , 
-oc JO 

where 

C ss (0 = (exp[*2n^(0]exp[-*2T^(0)])°ft(0 (B4) 

and 

C sQ (0 = C Qs (0=v 2 C d (0Gi(0 (B5) 
Cqq(0 = v*C dd (0Gi(0 , (B6) 
0,(0 = (^^(Oexp[i2n^(O]eaq>[-t2n0(O)])° 
C dd (0 = (a^(^9 x 0(O)exp[i2n ( /»(O]exp[- t 2n0(O)]) 



(£ is an abbreviation for (x,t)); 



(B7) 



In passing, we note that in general the presence of ap- 
proximate conservation laws implies the existence of large 
length scales like I2 on which transport is inhomogcncous. 
This may be related to the experimental observation 2 - 
that heat transport in spin chains is often extremely in- 
homogeneous. 



APPENDIX B: THE CALCULATIONS 



To evaluate the matrix elements Mi'Li, we first find ex- 
plicit expressions for the operators F* ml (t) = i[J p , Hnml\ 
(I = 0, 1), where H% m0 = H% m , H" ml = and the 

currents are given by Eqs. fZ5jl. We restrict our 

calculation to the even n terms (|12|l . (|19f) . The extension 
to odd n [Eqs. ((141) an d HOH is straightforward: due 
to the extra factors of d x <f>, it will essentially amount to 
trading C ss by C d d below [see Eqs. (JB4J), (|B6l) ]. We also 
define all operators below for a single chain. To leading 



s 

A QQ 

Eq. CEU) yields 



A, 



(2-Ka) 2n 

_ ±<,gV ml ?n*Kv 



{2tto) 



In 



4(g^,i) 2 " 2 
(27ra) 2 " 



G (x,i) = <^,t)0(O,O))° 



(B8) 



(B9) 



where at finite T the Green's function G<p(x,t) is given 



K 

G${x,t) = —In 



iraT/v 



sinh{7rT(a; — vt + ia)/v} 

naT/v 
sinh{7rT(ir + vt — ia)/v} 



(BIO) 
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Go(0 — lj an d the phonon propagator Gi(0 = G p (£) (at with B a numerical constant, 
finite T) has the form 



G p (x,t) « B 



iraT/vp 



sinh{7rT(x — v p t + ia)/v p } 
■naT I v p 



sinh{7rT(a; + v p t — ia)/v p } 



(Bll) 



J 



Substituting in Eq. (|B9I) and using the notation Bi, where Bq = 1 and B\ — B 1 we obtain 

2n 2 K / m\ 2/ 

C ss (x,t) =B l 



iraT \ ' 
v p ) 



sinh{7rT(a; - vt + ia) /v} sinh{7rT(a; + vt- ia) /v}} nK 

x [sinh{7rT(a; — v p t + ia)/v p } sinh{7rT(a; + v p t — ia)/v p }]~ 1 . (B12) 



Using the identity 



d x (j)(x,t) = lim lim(i7) d y exp[ijcf>(y, t)] 



7— >0 y— >x 

we can also express Cd, Cdd in terms of the function G^(x,t) and its derivatives: 

C d (x,t) = -2nd x G,p(x,t)C ss (x,t) 
C dd {x,t) = C ss {x,t)[An 2 {d x G^{x,t)) 2 +d 2 x G^{x,t)], 

where 

KttT, 



d x G^{x,t) 



4 v 



and 



■[coth{7rT(a; - vt + ia)/v} + coth{7rT(x + vt - ia)/v}] 



d x G$(x,t) = — — [sinh {ttT^ - ut + ia)/«} + sinh {7rT(a: + vt - ia)/v}}. 



K /ttT 



(B13) 
(B14) 

(B15) 
(B16) 



We now recall Eq. lj3Tjl for M^ ml , and consider the limit lu — > where the correlation functions F* given 



(B17) 



by Eq. (|B3D arc expanded to linear order in a;. We then get 

/OO /*00 
dxe- 4Ate / dttJm{C pq (x,t)}, 

where we have used the fact that C*(x,t) — C pq (x,—t) [see Eqs. I|B3|) through (|B16|) ] and hence the function 
tlm{C pq (x, t)} is symmetric with respect to t — > — i. The matrix element M*f(Afc, T) can be computed by substituting 
Eq. jEHl i a pl7| an d per fo rming the integrals. To compute M^f(Ak,T), M^ Q {Ak,T) we insert Eqs. (HHJ), l|B12)l . 
(|B15|) and (jB16|) into (|B13|) . l|BT4|l . This yields 



A sQ C s Q(:E,i) 



(27ra) 2 " a 



fsQ(x,t) 



(B18) 



f s Q(x,t) = [coth{TrT(x—vt+ia)/v}+coth{TrT(x-\-vt—ia)/v}] [sinh{7rT(a; — ui + ia)/u} sinh{7rT(x + — ia)/v}] 

x [sinh{7rT(x — v p t + ia) / v p } sinh{7rT(a; + v p t — ia) /v p }} 1 , 



-71 K 



and 



A QQ c QQ(x,t) 



(g^ ml ) 2 vWKB l fnaT\ 2n K+2 ( m T\ 
(2na) 2n a 2 



v p J 



fQQ(x,t), 



(B19) 



fQQ(x,t) = [n 2 ii'(coth{7rT(x - vt + ia)/v] + coth{7rT(a; + vt - ia)/v}) 2 + sinh 2 {ttT(x - vt + ia)/v} 
+ sinh~ 2 {7rT(.T + vt — ia)/v}~\ [sinh{7rT(a; — vt + ia)/v} sinh{7rT(a; + vt — ia)/v}} n K 

x [sinh{7rT(x — v p t + ia) / v p } smh{irT (x + v p t — ia)/v p }} 1 . 
I 

The last step is to insert the functions given by Eqs. the integrals. We employ the change of variables s = 
(|B12|) . llBT8)) and l|BT9)l into Eq. l|BT7|l and evaluate 
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T[x/v p — t), s' — T(x/v p +t), and define a dimensionless 
parameter A = v p \Ak\/T. Eq. (|B17|) is then recast in 
the form 

M^f (Afc, T) = T V Pq +2(n 2 K+l)-3 ( B2Q ) 

dsds'e- l ^ n( - Ak ^ s+s '^ 2 T^( S , s') 



oo poo 



— oo J — oo 



where r\ ss — 0, r\ S Q = 1, tjqq = 2 and the functions 
3~nl( s i s ') ^° n °t contain any dependence on Afc and T. 
Hence the double integral depends on them only through 
the parameters A and sgn(Afc) in the exponential factor. 
In particular, the high and low T limit cases are distin- 
guished by A -C 1 and A> 1, respectively. In the high T 
limit the exponential factor is expanded up to first order 
in A. The leading contribution to M^J and M®® comes 
from the O'th order, i.e. the integration results in a con- 
stant independent of Afc and T. However, since C s q(x, t) 
is an odd function of x, the leading contribution to M^f 
comes from the first order leading to an overall factor of 
Ak/T. These approximations yield Eq. (|42l) . 

To get the low T limit expressions for M^(Ak, T), we 
evaluate the integrals in Eq. (|B20(I in a saddle point ap- 
proximation where the large parameter is A. As implied 
by Eqs. (|BT2|) . ||BT8)| and (Efljty . the functions C pq (x,t) 
and hence !F^f(s,s') have branch-cut singularities along 
the imaginary axis (this is for the generic case where K 
is not an integer), one of which is close to the real axis. 
The integration over s' is performed first, with a slight 
deformation of the real axis to include the single saddle 
point 



A 



iO" 



(B21) 



Then, for sgn(Afc) < (sgn(Afc) > 0), the contour of 
integration over s is deformed from the real axis to the 
upper (lower) imaginary axis. The series of saddle points 
{si}°^ dominating this integral are close to the zeros of 
the sinh functions in (|B12|) up to a correction of order 
1/A — > 0. Since the contribution of each saddle point 
involves an exponential factor e _A ' Si '/ 2 , the overall result 
will be dominated by the minimal Si for which the func- 
tions JFP 9 (s,s') do not vanish by symmetry. The latter 
requirement excludes the contribution of sq = s' . The 
leading contribution therefore originates from the saddle 
point 



si = ±i . 



(B22) 



where sgn(Afc) = =p. Noting once again that v p <C v, we 
obtain the low T expression for M n i(Ak,T) in Eq. l|15fl. 
The matrix element M^(Ak,T) is then given directly 
up to a prefactor. To get the other matrix elements in 
(@2J>> we use scaling arguments noting that as long as 
(T/v\Ak\) <C (vp/v), the dominant momentum scale is 
Afc. 

The derivation of M«(Afc,T) [Eqs. through (|S7|l ] 
in the case where spinons and phonons couple linearly 



involves essentially the same calculation. In this case, 
however, the memory matrix is dominated by the contri- 
bution of the I = Umklapp terms Eq. {T^J, tH}. In 
addition, the spinon-phonon mixing introduces the fields 
4>i, (f>2 (the eigenmodes of the Hamiltonian H*), in terms 
of which the spinon field <f> can be written as 



>K | Cfa - S [ - 



1/2 



(B23) 



where we have used the definition = yK<p and the 
transformation Eq. (|53|l . As a consequence the operators 
F%m = i-[JpiHnm\ (f° r even n ) are given by 



Fnrn (fit 



(2na) n 
{2ira) n 



dx(e 



iAkx i2n<fi(x ,t) 



h.C.) 



dx(e iAkx e t2n4,{x ^ - h.c.) 



x v\Cd x 4>i - v 2 S 



1/2 



d 



xV2 



.(B24) 



The correlators (F% m ', F% m )^ are again written in terms 
of the functions C pq (x, t) [Eq. I|B3|) ]. Here, 

C ss (0 = (cxp[i2n0(O]cx P H2n0(O)])° (B25) 



and 

Cqq® 



\ 1/2 

v\ 1/2 

v$C'c u (0 - 2(v lV2 ycs [ — ) c 12 (o 



v„ 



(B26) 



{d x M0 eatp[i2n0(O] exp[-i2n0(O)])° 
(d x M09 x <pw (0) exp[*2n^(0] exp[-i2n^(0)]) 

and the coefficients A pq are given by 
4(gg n ) 2 (nKv) 2 



-4. 



A S Q = Aq s = 



(2na) 2n 

A{gV m ) 2 n 2 vK*/ 2 



(27TO) 



2n 



(B27) 



A QQ = 



±{g u nm ) 2 n 2 K 



(2ira) 2n 

Inserting Eq. [|B23j> into (|B25|) we get 



e 4n 2 KC 2 G!(x,t) e 4n 2 K S 2 (v/v p )G 2 {x,t) 



Css [x , t) 

G v {x,t) = (M^t)Mo,o)}° 



(B28) 



where similarly to Eq. IjBlOl 



G v (x,t) = -^l n 



iraT/v v 



4 



sinh{7rT(a; — v v t + ia)/v u } 
iraT/ v v 



'B29) 



sinh{7rT(a; + v v t — ia)/vv\ 
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Also, similarly to the derivation of Eqs. (|B13|) . (|B14|) we 
obtain expressions for C v , C vv i in terms of the functions 
G v {x,t) and their derivatives: 

Ci{x,t) = -2nVKCd x Gx(x,t)C ss {x,t) (B30) 

/ \ 1/2 

C 2 (x,t) = -2nVKS^—J d x G 2 {x,t)C ss (x,t), 



a n = n 2 KK 1 C' 1 and j3 n = n 2 KK 2 S 2 (v/v p ). Note that 
this expression is essentially the same function as (IB12|) 
with different parameters: v — > vi, v p — > v 2 , n 2 K — > a n 
and I — > (3 n . As a result, the integrals in Eq. I|B3(I can 
be evaluated in the same manner, yielding Eqs. I|64|) 
through (|6Tjl . 



C n {x,t) = C ss {x 1 t)[in 2 KC 2 {d x G 1 {x,t)) 2 +d 2 x G 1 {x,t)] 

1/2 



(B31) 



C 12 (x,t) = An 2 KCS — 



x d x Gi(x,t)d x G 2 (x,t)C ss (x,t) 
C 22 (x,t) = C ss (x,t) 

v 
u 



4n 2 KS 2 — (d x G 2 (x, t)f + d 2 G 2 {x, t) 



The explicit dependence on x and t is obtained from Eqs. 
(|B29|) . l|B28)l . which yield 



[sinh{7rT(x — Vit + ia) /v{\ sinh{7rT(a; + V\t — ia) /i>i}] a " 
[smh{nT(x—V2t+ia) /v 2 } sinh{7rT (x+v 2 t—ia) /v 2 }} 

(B32) 



APPENDIX C: ERRATUM 



During the preparation of a subsequent article^ we have spotted an error in the derivation of the effective propa- 
gator of 3-dimensional phonons G p (t, x), obtained after integration over the momenta in the direction perpendicular 
to the spin chains. The revised calculation of the thermal conductivity n(h — 0) as a function of T [Eq. (47)] yields a 
modified power-law prefactor, however, our main result - the exponential behavior with the characteristic temperature 
scale T* [Eq. (48)] - is unchanged. 

The effective single-phonon propagator G p (x,t) = (d x q(x,t)d x q(0,0)} (where q is the displacement field in a par- 
ticular chain) can be derived by an inverse Fourier transformation of G p (k,t) — k 2 J d 2 k±G p D (k,t), where k is the 
component of k along the chain and G p D (k,t) the correlator (<7k(*)<?k(0)) of free (isotropic) 3-D-phonons at finite T: 



G 3D (k,t > 0) = const x 



kit 



+ (1+rikK 



kit 



(CI) 



where rtk = (e" p ' k '/ T 
Appendix B: 



1) 1 is the phonon occupation. We obtain the following expression, to replace Eq. (Bll) in 



G p (x,t)=B[ 



± ^ * 

v p t \ {x + v p t) 3 (x — v p t) 3 

f) EE 



(n + critT)(n + ai(t — x/v p )T) 3 [n + aitT)(n + ai(t + x/v p )T) 3 



(C2) 



with B a numerical constant. 

To find the contribution M^f(Ak,T) to the memory matrix (where / ^ denotes the number of phonons involved 
in the Umklapp process), we use Gi{£,) = [G p (£,)] 1 in Eq. (B9) with G p given by Eq. (H21> above. At low T, the 
single-phonon process / = 1 yields a subdominant contribution ~ e - v P \ Ak \/ T _ The physical reason for that is that the 
energy cost of the process involves an excitation of a single phonon that carries the entire momentum transfer Afc, 
i.e. v p I Afc I . It is energetically favorable to distribute the momentum equally between the initial and final state using 
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a two-phonon process (I = 2): a phonon with momentum Ak/2 is scattered to a final state with momentum — Ak/2. 
This yields the leading contribution to the matrix clement M TT : 



M[ 2 T ~T 2K + 3 exp 



v p \Ak\ 



2T 



(C3) 



which should replace M^ T in Eq. (46) (with Afc = G/2). Higher order processes (I > 2) give the same exponential 
dependence but are supressed by power law prefactors due to phase space restrictions. The resulting expression for 
the thermal conductivity is 



n(h = 0) ~ k,q 



2K 




J exp 


T 



~ 4 



(C4) 



This should replace Eq. (47). 
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